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Abstract 

In commodity markets, the convergence of futures towards spot prices as the time to 
maturity of the contracts goes to zero is usually justified by no-arbitrage arguments. In 
this paper we propose an alternative approach, that relies on the expected profit maxi¬ 
mization problem of an agent producing and storing a commodity while trading in the 
associated futures contracts. In this framework, the relation between the spot and the 
futures prices holds through the well-posedness of the maximization problem. We show 
that the futures price can still be seen as the risk-neutral expectation of the spot price 
at maturity and we propose an explicit formula for the forward volatility. Moreover, 
we provide an heuristic analysis of the optimal solution for the production / storage 
/ trading problem, in a Markovian setting. This approach is particularly interesting 
in the case of energy commodity: it remains suitable for commodities characterized by 
storability constraints, when standard no-arbitrage arguments can not be safely applied. 
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1 Introduction 


In this article we aim at explaining, through a parsimonious model, the relation between the 
spot and the futures prices of a commodity, that can be storable or not. The prices relation 
derives from the well-posedness of the simple optimization problem of an operator involved 
in the production and (when possible) the storage of a commodity. This producer also trade 
in the futures market of the commodity. 

Such an approach is different from the classical no-arbitrage reasoning usually employed 
to explain the price relations in commodity markets. Interestingly, it remains relevant for 
commodities characterized by storability constraints, when standard no-arbitrage arguments 
can not be safely applied. This does not mean, however, that our framework is devoted 
only to non-storable commodities, like electricity, nor that it is focused on commodities with 
storable inputs (Aid et al , 2009 and 2013 [1, 3]). 

Commodity futures prices fail to satisfy the no-arbitrage relation that holds for investment 
assets such as bonds and stocks, which can be stored at no warehousing or depreciation costs 
(Eydeland, 2002 [23], Chap. 4, pp. 140-43). This no-arbitrage relation states that the futures 
price F(t,T ) of a contract written on a stock, observed at date t for a delivery at T, is the 
spot price St capitalized at the interest rate r between t and T: F(t,T ) = Ste r ^ I '~ t \ 
Introducing a storage cost in the analysis is however not sufficient to depict the behavior 
of the prices spread in commodity markets. There is still a need to explain, first negative 
spreads 1 , second the simultaneous presence of positive inventories and backwardated prices 
(see among others, Working 1933 [41]). The no-arbitrage relation holding for stocks and 
bonds can only be preserved by the introduction of an extra variable: the convenience yield 
(see Kaldor, 1940 [26] and Lautier, 2009 [29]). The latter represents an implicit revenue. 
The presence of such a yield, which is associated to the physical commodity but not to the 
futures contract, explains that the operators maintain their inventories even in the presence 
of negative spreads. With the introduction of such a variable, the former relation becomes 
F(t,T ) = Ste^ r ~ y ^ T ~ t \ where y is the convenience yield net of the storage costs. 

An important effort has been undertaken, in the economic and financial literature, in order to 
propose models that rely on storage behavior to explain and reflect the interaction between 
the spot and futures prices of a commodity. Within this abundant literature, our article 
is close to those of Brennan (1958 [14]) and Routledge et al (2000 [37]). In these works, 
the authors develop production-storage models that connect the spot and the futures prices. 
Nevertheless, the models are mainly developed in order to allow for comparative statics and 
little information on the conditions under which no-arbitrage holds is given (this is especially 
true for the model of Routledge et al (2000 [37]). 

Apart from these models, pure financial models explicitly designed for the pricing of com¬ 
modity derivatives and relying on the concept of convenience yield were also built (Lautier 
(2005 [28]). In our opinion, the most important reference, in this literature, is probably 
Schwartz (1997 [39]). In this setting however, the relation between the spot and the futures 
prices relies on the hypothesis that there exists a unique risk-neutral measure that one can 
calibrate thanks to market data. 

The development of electricity markets in the last thirty years has challenged the idea that 
the futures price of a commodity is linked to its spot price by a convenience yield. Electricity 
indeed cannot be stored. Still, futures prices exhibit both contango and backwardation struc¬ 
tures (Benth et al, 2013 [11]). Many authors have pointed out that, in the case of electricity, 
the convenience yield may not apply. The same is true for the no-arbitrage method used 
in mathematical finance to obtain a risk-neutral measure (Benth et al, 2003 [7]). Moreover, 

X A positive spread (a contango) arises when at date t, the futures price F{t,T) is higher than the spot 
price St- A negative spread (a backwardation) is a situation where F(t,T ) < St- 
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due to the particular nature of the futures contracts negociated on these markets (which are 
basically swaps, see Frestad, 2010 [10]), even the convergence of the futures price to the spot 
price is an issue (Viehmann, 2011 [40]). 

Since the pricing of derivatives can not rely on the concept of a convenience yield in the 
case of electricity, this non-storable commodity has fostered researches on how to restore a 
relation between the spot and the futures prices. The first approaches relied on two-date 
equilibrium models (Anderson and Hu, 2008 [5], Bessembinder and Lemmon, 2002 [12], Aid 
et al, 2011 [2]), where the risk-neutral measure is extracted from the risk aversion parame¬ 
ters of the agents. A more complex approach consists in the extension of the market beyond 
the underlying asset, to include production factors (fuels) (Aid et al, 2009 and 2013 [1, 3]), 
production constraints (Bouchard and Nguyen Huu, 2013 [13]) or gas storage levels (Douglas 
and Popova, 2008 [22]). In this setting, the spot price is the result of an equilibrium between 
production and consumption. The futures price is obtained as an expectation of the spot 
price. 

In this article, we perform an analysis of how the futures prices can be related to the spot 
prices by arbitrage arguments, independently of the storability properties of the underlying 
asset. We do not intend, however, to provide a model especially designed for non-storable 
commodities. This kind of work is performed in Benth et al (2007 [8] and 2008 [9]), Meyers- 
Brandis and Tankov (2008 [31]), and Hess (2013 [24]), among others. Nor do we aim at 
proposing a general theory of the joint dynamics of the spot and futures prices. 
Consequently, this framework could be applied to any commodity, whether it is storable or 
not. This could be useful for a large range of industrial companies operating in commodity 
futures markets, such as energy utilities, airline companies, processors of steel or of agricul¬ 
tural products. The profit maximization of such companies indeed relies on the spot prices 
as well as on the futures prices. Within our framework, once a spot price model and the 
market price of demand risk are chosen, there is only one futures price model that would be 
consistent with their optimization process. 

The remaining of this article is organized as follows. In Section 2, we explain how we obtain 
the no-arbitrage relation through a simple model of production, storage and trading, where 
an operator maximizes his expected utility. This operator has no impact on the spot price. 
The storage cost and the production functions are supposed to be convex. The storage ca¬ 
pacity is bounded as well as the instantaneous storage or withdrawal. Moreover, the operator 
has access to a derivative market where a futures written on the commodity is negotiated. 
The contract is traded on one maturity only and the liquidity of the futures markets is sup¬ 
posed to be unlimited. Naturally, the operator has no impact on the futures price. 

In Section 3, inspired by previous results established in, e.g., Rogers (1994, [36]) and Ankirch- 
ner and Impeller (2005 [6]) for a pure trader, we show that the existence of a risk-neutral 
measure is the consequence of the finiteness of the operator’s value function. Stated dif¬ 
ferently, this result means that if there were no risk-neutral measures, the operator could 
take advantage of his production capacity or storage facilities to get an infinite utility. In 
the same section, we also prove that the futures price always converges to the spot price, 
regardless of the storability properties of the commodity. This result is important, especially 
for the analysis of electricity markets. It has been indeed pointed out that the electricity 
futures prices predict realized spot prices rather poorly (Prevot et al, 2004 [34]). This study 
was however done on monthly contracts while the convergence issue concerns only maturities 
close to zero. For instance, day-ahead futures contracts quoted on the German electricity 
market exhibit lower discrepancy with the realized spot prices (see Viehmann, 2011 [40]). 
We finally discuss in Section 4 the trading-production problem faced by the operator, with 
a specification of the demand dynamics. We obtain an explicit formula for the volatility of 
the futures contract and we relate it to the volatility of the underlying conditional demand 
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for the commodity. Moreover we argue that, in a Markovian setting and for an agent having 
a power type utility function, the optimal command for the management of storage is of a 
bang-bang type. Not surprisingly, the decision to store or to withdraw the commodity is 
based on the comparison between the spot price and the ratio between the marginal utility 
of one unit of storage and the marginal utility of the wealth of the operator. Lastly, in this 
setting the optimal trading strategy on the futures market is such that the operator holds a 
long position when the futures prices exhibit a positive trend. 


2 The model for the individual producer 

In this section we explain how we obtain the no-arbitrage relation through a simple model 
of production, storage and trading. 

Let (O, P) be a filtered probability space satisfying the usual conditions, i.e. (fFt) 

is a P-completed and right-continuous filtration. All processes considered in this paper are 
assumed to be defined in this space and adapted to this filtration. 

Our agent is the producer of a commodity who also acts in the derivative market associated 
to this commodity. This producer is a price taker. On the physical market, he has the 
possibility to sell all his production, to store a part of it or to reduce his stocks. On the 
derivative market, he has the possibility to buy or to sell a certain amount of a unique futures 
contract. These choices depend on the price conditions he faces, both on the physical and 
on the derivative markets. In what follows, we will first focus on the physical market. Then 
we will expose the trading activity on the derivative market. 

2.1 The profit on the physical market 

On the physical market, the agent has the possibility to decide how much he produces and 
to manage his stocks dynamically. His instantaneous profit TTt can be written as follows: 

7T t = (q t - u t )S t - c(q t ) - k(X t ) (2.1) 


where: 

• qt is the production of the agent, 

• Ut is the amount stored ( ut > 0) or withdrawn ( ut < 0), 

• St is the spot price of the commodity, 

• c : M+ —)• R is the production function, 

• k : M+ —)• M is the storage function, with k( 0) = 0, 

• Xt is the storage level at time t. 

Assumption 2.1. We assume that both functions c and k are differentiable, strictly increas¬ 
ing, strictly convex and nonnegative. 

Moreover, we will always work under this standing assumption on the spot dynamics: 
Assumption 2.2. Let (St) be a bounded continuous process. 

The constraints faced by the producer on the physical market can be summarized as follows: 
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• the agent’s production cannot exceed his capacity q : qt £ [0, q] for some q > 0; 

• instantaneous storage and withdrawal are bounded, with ut £ \u, u] for given thresholds 
u<0<u] 

• the storage capacity itself is bounded by some X > 0, so that adding the positivity 
constraint on the inventories we have Xt £ [0, X] a.s. for all t £ [0,T]; 

• the storage dynamics is: dXt = Utdt, with Xq = uq > 0. 

In this setting, there is no uncertainty on the production of the commodity. The only source 
of uncertainty comes from the demand side. This hypothesis is made in order to simplify the 
computation and the formulas. The introduction of a stochastic production capacity driven 
by a Brownian motion independent of the one driving the demand would introduce market 
incompleteness. Besides, this hypothesis is reasonable for a large number of commodities, 
especially energy commodities. 

2.2 The trading activity on the derivative market 

We assume, for the sake of simplicity, that the interest rate is zero. Moreover, there is only 
one futures contract available, for a given maturity T > 0. The price of this contract at t is 
Ft = Ft(T). We assume that the futures price process (Ft) is a continuous semi-martingale, 
adapted to the filtration (Ft). The trading portfolio on this contract is given by: 

r T 

V$= / 9 t dF t , (2.2) 

Jo 

where 8 is a real-valued predictable (T))-integrable process. We assume that the futures 
market is liquid, which is standard in this context. On the empirical point of view, this is 
the case for the short-term maturities of commodities like crude oil and natural gas. This 
is however less obvious, at least up to now, for electricity markets (A'id, 2015 [4], Chap. 2, 
Sect. 2.2.3). 

2.3 The production-trading problem 

As a commodity producer, the agent acts so as to maximize the expected utility of his 
terminal wealth. His utility function U : —>• [—oo,oo[ satisfies Inada conditions and 

is such that U(x) —> oo whenever x —>• oo. Moreover, we assume Reasonable Asymptotic 
Elasticity (RAE: see Kramkov and Schachermayer, 1999 [27]): 

AE(U) := lim sup < L 

x^oo U(x) 

In this setting, we propose the following production-trading problem: 

r(r 0 ):=supE u(tq+ [ TT t dt + V? + 8t~(Ft — St)) , (2.3) 

u,q,6 \ Jo J _ 

where: 

• ro > 0 is the initial wealth of the operator, 
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• 7 Tt is the instantaneous profit on the physical market, expressed as a function of the 
quantities produced and stored, given by equation (2.1), 

• Vf is the trading portfolio in the futures market as given by equation (2.2), 

• the term Ot~(Ft — St) can be explained by the delivery conditions of the futures 
contract at expiration (see the heuristic discussion in discrete time below). 

The controls (it, q, 0) have to satisfy the following additional constraints: 

• constraint on the wealth of the agent to prevent infinite borrowing: 

R X t ' U ’ q ’ 6 ■= r 0 + f Ksds + V t e + 9 T -{F t - S T )l {t =T) > 0, t G [0, T\. (2.4) 

Jo 

• the production-storage controls (ut, qt)te[o,T] are adapted processes with respect to the 
filtration (Ft) and they satisfy the constraints previously described in Paragraph 2.1. 

Discrete-time heuristics. We provide some heuristics in discrete-time to justify the form 
of the continuous-time problem (2.3). The terminal total wealth for a trader-producer who 
produces, e.g., energy out of fuels, and trades in futures contracts on energy over the finite 
time grid {0,1,... , T} with T € N can be written as follows: 

T—l T- 1 

V T = J2l(qt- u t) s t- c (qt)-k(X t )\ + ^ 9 t (F t+1 - F t ) 

t= 0 t=0 

+ 6t—\Ft—i — hr St — c(qx) — k(Xx), 

where hr is the quantity bought or sold at terminal date that allows to fulfill the commitment 
taken on the futures market at the expiration of the contract, i.e. h t is such that 

0t —i = hx + qr — ut, 

so that the terminal total wealth becomes 

T T—2 

Vt = 53(9* " u t) S t ~ <qt) ~ k(X t )\ + °t{F t+ 1 - F t ) + 0 T -i(F t - S T ), (2.5) 

t =0 t =0 

which constitutes the discrete-time analogue of the total wealth appearing inside the utility 
function in the equation (2.3). 

In the next section, we will prove that a necessary condition for the problem (2.3) to be 
well-posed is the equality Ft = St, so that the third summand in the equation (2.5) (as well 
as its continuous-time counterpart in (2.3)) vanishes. 

Our first objective is to deduce, from the well-posedness of the problem described by equation 
(2.3), the no-arbitrage condition that would link spot and futures prices, i.e. to prove that 
there exists some equivalent probability measure Q such that 

F t = E Q [S T \F t \, *G[0,T], 
and to try to compute explicitly this futures price. 
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3 Existence of the optimum and spot-futures no-arbitrage re¬ 
lations 

In this section we derive the existence and the uniqueness of an optimal solution ( q*,u*,9*) 
for the optimization problem (2.3). At the same time, we obtain no-arbitrage relations 
between the spot and futures prices, as well as the convergence of the futures prices to the 
spot prices when the time-to-maturity goes to zero. Moreover, we show that the optimal 
production q* can be computed explicitly even in this quite general framework. The way 
we obtain an explicit expression for the other optimal quantities, i.e. storage and trading 
activity ( u*,6 *), is discussed in Section 4. 

Assumption 3.1. Let v(ro) < oo for some initial wealth ro > 0. 

Convergence of the futures to the spot prices and no-arbitrage relation. Our 

first result states that, as long as our optimization problem is well-posed, one must have a 
convergence of the futures price towards the spot price when the time-to-maturity tends to 
zero. This is true even when the underlying commodity of the contract is non storable. 

Proposition 3.2. Under Assumption 3.1, we have Ft = St- 

Proof. Assume that P (Ft ^ St) > 0 and let A = {Ft > St} and B = {Ft < St}- Consider 
the following trading-production strategy: 

q = u = 0, 6 t := (al A ~ fdl B )lt=T, 

where a and (3 are arbitrary positive numbers. Since A and B are .^--measurable (St and Ft 
are both continuous processes), 8 is clearly a predictable and F)-integrable trading strategy. 
Pursuing such a strategy gives terminal wealth: 

x + oi(Ft — St)1a + P(St — Ft)1b, 

so that, letting a -A- oo and fd = 0, if P(A) > 0 or /? —> oo and a = 0 if P(B) > 0 we get 
v(x) = oo (recall that U(x ) —>• oo when x -A oo), which contradicts the well-posedness of our 
maximization problem. Thus, we can conclude that a.s. Ft = St- □ 

Before proving the existence of a solution to our optimization problem, we deduce the no¬ 
arbitrage property for the futures contracts from the finiteness of u(ro). This is the content 
of the next proposition, which adapts arguments from Proposition 1.2 in Ankirchner and 
Imkeller (2005 [6]), where No Free Lunch with Vanishing Risk (henceforth NFLVR) is deduced 
from the well-posedness of an optimal pure investment problem. We refer to Delbaen and 
Schachermayer’s paper (1994, [20]) for its definition as well as the proof of their celebrated 
version of the fundamental theorem of asset pricing. In this paper, we deduce something 
stronger, namely a variant of NFLVR not only for the trading portfolio but also for production 
and storage. To be more precise, let us redefine the NFLVR condition for our setting. 

Definition 3.3. A Free Lunch with Vanishing Risk with production and storage is a sequence 
of admissible plans (qf,uf,9f), n > 1, such that R^f n := R^f q ,u converges a.s. towards 
some nonnegative r.v. Rj- satisfying P(i?^ > 0) > 0 and ||(i?^,) - || 0O —> 0 as n —>• oo. We 
will say that NFLVR with production and storage is satisfied if there are no such strategies 
in the model. 

Notice that since the production and storage controls are bounded, there exists a constant 
M > 0 such that | 7Tt.dt\ < M for any admissible (q,u) giving the instantaneous profit 7q. 
This fact will be used in the proof of the following result. 
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Proposition 3.4. Under Assumptions 2.1, 2.2, 3.1 and for all vq > M we have that v(ro) < 
oo implies NFLVR with production and storage. In particular, NFLVR for futures prices (Ft) 
also holds, which in turn implies that there exists at least one equivalent probability measure 
Q under which Ft is a local m.artingale. 


Proof. Suppose that NFLVR with production and storage is violated, so that we can find 
a sequence of terminal payoffs R'f = irfdt + J Q r OfdFt such that R'f —>• Rf for some 
nonnegative r.v. Rf with P(i?§, > 0) > 0 and e n := \\(R'f)~\\oo —>• 0 as n —>• oo. Here ir n 
denotes the instantaneous profit coming from a production q n and a storage u n . Now, take 
any ro > 0. There exists a 6 > 0 such that k := U(vq — M — 6) A 0 > —oo. Consider the 
sequence of strategies ( q n ,u n ,9 n ) where 

q n = — q n /\q, u n = — /AmVH), 6 n = — 9 n , 

and let R n denote the corresponding payoff. Clearly, R'f > — 6 and moreover the r.v. C/(?’o + 
Rf) is bounded from below by k. Indeed: 

U fro + £ irfdt + Vf \ = ufr 0 + £ nfdt + ^ v t \ > U (r 0 -M + = « > 


Since Rif converges to the nontrivial nonnegative r.v. Rf and any profit is bounded by a 
constant M, one can find an integer no and real numbers b,c > 0 such that P(V^T > b) > c 
for all n > no- Since k < 0 we have 


lim inf E 


U ( r 0 + / 7 tfdt + Vf 

Jo 


= lim inf E 

n 

> lim inf E 

Tl 

> lim inf E 


U\r Q + \ T nr?dt+—Vf 

0 £ n 

■T 

kS 


A{v™<b} + U ^ro + J n'fdt + —6^ 1 


k( 1 — c) + U [ vq — M - b ) c 

£yi 


= OO, 


yielding that v(r o) = oo. To conclude, just recall that NVLFR with production and storage 
clearly implies NFLVR for trading in futures contracts, which yields in turn the existence 
of at least one equivalent probability measure Q for which Ft is a local martingale (see the 
seminal paper of Delbaen and Schachermayer, 1994 [20]). □ 


Remark 3.5. The previous result implies in particular that Ft is a local Q-martingale under 
some equivalent probability measure Q with terminal value at time T given by the spot price 
St- If Ft was a true martingale under Q, we would recover the well-known relation 

F t = E q [S t | F t ], t € [0, T], (3.1) 

This would be the case if, e.g., Ft was an Ito process whose Brownian part had a square- 
integrable volatility process (as in the next section), or more generally if Ft was of class (D) 
(see, e.g. Protter 2004 [35, III.3]). 


Existence and separation principle. An immediate consequence of the previous con¬ 
vergence result is the following separation principle, stating that solving our optimization 
problem is equivalent to maximize first with respect to the production control q and then 
with respect to the storage and trading controls (■ u,0 ). On the other hand, maximizing the 
production can also be performed in two steps. Let us denote 


v(r 0 ) 


supE 

u,q,9 


u (V 0 + yf + zf + Vr'j 
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—oo. 


{ Vp>b} 



where we set 


Yf := [ T (q t S t - c(q t ))dt, Zf := - f \u t S t + k(X t ))dt. 

Jo Jo 

We can solve our problem in two separate steps. First we solve u(?’o) with respect to the 
production control q (for given it, 9 ); second, we solve with respect to the controls (u , 9). Let 
us start from the production side. 

Proposition 3.6. Under Assumptions 2.1, 2.2, 3.1, for any given admissible investment 
strategy 9 and storage policy u, the optimal production control q* is given by 

Qt = ( c/ ) _1 ('S't) 1 {(c')- 1 (St)<g}> * G [°’ T ]- ( 3 - 2 ) 

Proof. It suffices to maximize w-wise inside the integral in the term Yf containing the pro¬ 
duction controls. Differentiating with respect to (p for a fixed t gives St — c'(qt) = 0 so that, 
taking into account the constraint qt G [0 ,q] and since c is strictly convex, we have (3.2). 
The proof is completed. □ 


Let us denote 


Yf := Yff = / (qt.S t — c(q*))dt 


where q f is given by (3.2). Now, let us consider the optimal storage/trading problem 


v(r o) := supE 

u,6 


U ( To + Yj> + Zj' + Vrp 


(3.3) 


The next result establishes existence of a unique optimal storage/trading policy (u*,9*). 

Proposition 3.7. Under Assumptions 2.1, 2.2 and 3.1, there exists a unique solution (u *, 9*) 
to the problem (3.3). 

Proof. First of all, if one admits the existence of a solution, its uniqueness follows at once 
from the strict concavity of the utility function U. Let ( u n ,9 n ) be a maximizing admissible 
sequence for the problem (3.3), i.e. E [U (?’o + Yf + Zf + Vf)] —>• u(?’o) as ra —oo, where we 
denoted 

Zf := - [ T (u?S t + k(Xf))dt , X? := u 0 + 

Jo 

We prove the compactness property of the sequences u n and 9 n separately. 

For the sequence of storage strategies u n , we use the Komlos theorem, stating that for any 
sequence of r.v.’s (f n ) bounded in L 1 , one can extract a subsequence (£ nfc ) converging a.s. in 
Cesaro sense to a random variable € L 1 (see, e.g., Theorem 5.2 in Kabanov and Safarian, 
2009 [25]). We apply this theorem to the sequence of processes u n , that can be viewed as 
random variables defined on the product space (fl x [0, T], V, dPdt) where V is the predictable 
e-field. The sequence u n is clearly in L 1 since it takes values in the interval [—u, u]. Thus, 
there exists a predictable process u° taking values in the same interval, such that the Cesaro 
mean sequence u n := (1/n) converges a.e. towards u°. Indeed it is immediate to 

check that the sequence u n takes values in [—' u,u] as well. Moreover, the cumulated storage 
process along the new sequence, Xf := fg vJfds, is well-defined since each u n is bounded and 
it takes values in [0, X], By Lebesgue dominated convergence we have Xf —> X® := fg u®ds 


u n s ds , 


VP := 


9 n dFf. 
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a.s. for all t £ [0,T]. Since the function k is continuous, we have k(X™) —>• k(X °) a.s. for all 
t. Finally, thanks once more to the boundedness of the controls and to the continuity of k, we 
have \Zj,\ < C\ Sf/lt], which is bounded (since St is bounded uniformly in t). Therefore, 
applying the dominated convergence theorem again, we get —>• Zj, a.s. as n —> oo. 

As for the compactness of the sequence of trading strategies 9 n , we can work with the 
corresponding wealth process Cesaro mean sequence, that we denote by VTp. The admissibility 
property and the uniform boundedness of Z'j, yields that this sequence is uniformly bounded 
from below by some constant. Therefore, we can apply Theorem 15.4.10 in [21], implying that 
there exists a convex combination Vp £ conv(Vp, Vp +1 , ■ ■ •), which converges a.s. and whose 
limit is dominated by some Vp := O^dFf for some admissible 9°. Moreover, applying 
this procedure to the Cesaro means of storage strategies u n one gets another sequence of 
admissible storage strategies u n converging a.s. to the same process n° as u n before. 

To conclude the proof, we need to show that 

v(r 0 ) < E [U(x + Yp + V$ + Z%)\. 

To do so, it suffices to use the assumption that U satisfies RAE by proceeding as in the proof 
of, e.g., Theorem 7.3.4 in Pham (2000 [33]). Repeating his arguments gives us the inequality 
above getting that (u°,9°) is the optimal storage control (u*,Q*). The proof of existence is 
now completed. □ 

Remark 3.8. The fact that the trading-production problem above can be solved in successive 
steps does not mean that the optimal controls are independent. Only the production control 
q can be deduced independently from u and 9. Indeed, since the producer has no impact 
on the spot price, his optimal strategy is simply to equal his marginal cost of production 
with the spot price. Thus, no matter whether there exists a futures market or not, one 
would observe the same production level q. This is not the case for the optimal storage 
policy u and the optimal trading strategy 9: they are not independent. This fact has two 
consequences. First, the introduction of a futures market modifies the way storage capacities 
are managed. This point may be of interest for the econometric analysis of the relations 
between the level of the inventories and the prices. Second, as soon as the industrial process 
includes storage activities, the trading can not be separated from the storage without a loss 
of value. This should be taken into account for the organization of the trading activities in 
industrial companies. 


4 The optimal production-trading problem 

In this section we illustrate, in a simple setting, how the approach developed in the previous 
sections can lead to a consistent model for futures prices. 

We discuss the trading-production problem faced by the operator with a specification of the 
demand dynamics. Then we determine the volatility of the futures contract and we relate it 
to the volatility of the underlying conditional demand. 

4.1 The dynamics of the demand for the commodity 

In order to illustrate the previous approach, we use a simple dynamics for the demand. In 
particular, since our goal is not to propose a model that would match all stylized facts, we 
do not include features such as seasonality and jumps. 
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In commodity markets, the spot price St results from the availability of the raw material on 
the physical market. In our model, this availability is measured through the confrontation of 
the total capacities of the market and the demand for the commodity, in the following way: 

S t = b ■ g(C — D t ) ■ f (D t ), (4.1) 

with b a constant of normalization for dimension purposes, C > 0 the maximum available 
production and storage capacities of the market (supposed constant), Dt the total exogenous 
demand for the commodity (which is an (J^-adapted continuous process) and f(D ) the 
marginal cost of production and storage for a demand level D. 

As there is a non negativity constraint on inventories, the spot price can jump to very high 
levels when the total capacities are not sufficient to fully satisfy the demand. This behavior 
is captured by the scarcity function g: 

g(x) = l x>0 • min(l/ x, 1/e) + l x<0 \/e. 


The effect of scarcity on commodity prices is clearly illustrated for the case of oil in Buyuk§ahin 
et al (2008 [15] p. 56, fig. 10). The specific form of g above has been successfully implemented 
in the case of electricity spot prices in Aid et al (2013 [3]). 


Since the production optimization problem has been solved in Proposition (3.6), it remains 
to treat: 


v(x) = supE U (x + Yff + z t + V$) 

ufi L ' ' 


where 

Yt= [ (Qt S t ~ c(q*))dt, te[0,T}. 

Jo 

with q* as in (3.2). We recall that Zj, = uq + ut.dt is the cumulated storage and that 
Vrf- = / T 0 t dF t is the portfolio traded over the period [0, T], 

We assume that the futures price process (Ft) is an Ito process (or a diffusion to exploit the 
Markov property). More precisely: 


Assumption 4.1. 1. Let the demand for energy Dt be mean reverting, with a long-run 

mean set to zero: 


dD t = aD t dt + adWt, (4.2) 

where a, a are constants and W is a standard Brownian m.otion. We denote by (Ft) 
the natural filtration generated by it and completed with the P -null sets. 

2. Assume that the futures price F is an ltd process fulfilling 


dF t = a t dt + /3 t dW t , 


where ct,f3 are some (Ft)-predictable real-valued processes such that a.s. 



at\dt + E 



< oo. 


The integrability assumption on the volatility is here only to have the (true) martingale 
property of Ft and consequently the very useful formula Ft = E^ [St \ Ft ]. 


11 



4.2 Equivalent martingale measures and forward volatility 

First of all, we notice that Remark 3.5 and Assumption 4.1 imply that 

F t = E q [S t | F t \, t G [0, T], 

where Q is an equivalent martingale measure for the futures process (Ft), which means that 
Q must satisfy: 

L ‘ '= ^ l*= exp {- i x - dW - -1 i ■ 

where A is a (J-))-adapted process (viewed as “market price of demand risk”) such that 

• at — A tfit = 0 a.e. dF <g) dt, 

• Jo X l ds < °°> 

• E[Lj.] = 1. 

At this point, in order to specify completely the dynamics of the futures price under P, we 
need to assume a particular and tractable form for the market price of demand risk A t.. 

Assumption 4.2. Let us assume that Aj = \o(t) + Xi(t)Dt, t G [0, T], where Ao,Ai : [0, T] —>• 

M are deterministic functions such that the last three properties above are satisfied. 

A consequence of this assumption is that the drift at of the futures price takes the form 
a t = (Ao (t) + Ai {fi)Dt)fiti which is completely determined up to the volatility fit- 

We will see in what follows that the special form of the production function defining the spot 
price St in (4.1) implies a particular functional form for the volatility of the futures price 
process (Ft). 

Let Q be the equivalent martingale measure corresponding to the market price of demand risk 
A t as in Assumption 4.2. Under such a measure, the demand has a dynamics characterized 
as follows: 

dD t = ((a + Xi(t)a)D t + A 0 (t)a)dt + adW ^, 

where W (X is a standard Q-BM. Thus, the conditional distribution of Dt given Dt under Q 
is Gaussian with conditional mean mf T and variance given by 

mf T = etf (°+ Al ^ a)ds ^ Dt + e~ fo( a +^( u )^ du A 0 (s)a d s^ , = u 2 J^ e - 2 f t fi a +\i(u)<r)du ds 

(4.3) 

To complete the description of our model, we set a specific shape for the marginal cost of 
production. 

Assumption 4.3. Let the marginal cost of production f be equal to 

f(d) = d a l( 0 < d < M ' ) +M a l( d > M y d G M, 

for some exponent a £ (0,1) and some upper bound M > 0 such that our conditions on f 
are fulfilled. Moreover, let M > C — e. 
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Under all these assumptions, we can express the spot price St as a function of the demand 
St = if{Dt), where the function is given as follows: 


/ d a d a M \ 

V’(d) = b • ( — l(0<d<C-e) + (J __ d \c-e<d<M) + q _ ^(d>M) J • (4-4) 

Notice that the spot price St is always nonnegative. 

The futures price at time t computed under the above measure Q is given by 

F t = E?[tf(D T )}, t € [0,T], 

where Ef denotes the conditional Q-expectation given Ft = F t D . We denote by hr D t {y) the 
conditional density of Dt given Dt, that is: 


h'T,D t {y) 


i ( (>/ - ) j \ 

E t , T ^ eX U 2E It ) 


where the mean m9 T and the variance are given in (4.3). We recall that the variance 
does not depend on Dt- 

We can express the futures price Ft as a function of the demand at time t, Dt, as: 


F t = <p(t,D t ) 


i>{y)h T ,D t {y)dy. 


A simple application of Ito’s formula together with the martingale property of the futures 
price Ft under Q gives that the volatility of the futures price, f3(t, Dt), is given by: 

Pt = ff = a^(t,D t ). 

If we compute explicitly the first and second derivatives of the futures price, c p(t,Df ), with 
respect to the demand, we obtain the following result giving a complete specification of the 
parameters of the forward dynamics. 

Proposition 4.4. Under Assumptions 4-1, 4-® an d 4-3, the well-posedness of the optimal 
production-trading problem (2.3) implies that 

dF t = a t dt + P t dW t 


where 


a t = a(t, D t ) = (A 0 (t) + Ai (t)D t )/3 t , 

Q 

Pt = P(t, D t ) = a I ^(y) V ~™* |T ef< T ^ +x ^) du h T)Dt (y)dy, 

Jr *- l t,T 


for all t € [0, T]. Moreover, the forward volatility (3{t,Dt) is increasing in the demand. 


4.3 The production-trading optimization problem in a Markovian setting 

In this section, we provide for the sake of completeness an informal discussion of the optimal 
solutions within the Markovian model determined in the previous proposition. Let us assume 
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that the preferences of the agent are of power type, i.e. U(x) = x 1 , x > 0, where 7 G (0,1). 
Recall that the problem we want to solve is the following: 


v(x) 


sup E 
( u,q,6)eA 



7 Tfdt + 


(4.5) 


where ?’o > 0 is the initial wealth. Recall that 77 is the profit rate given by: 


77 = ((It - u t )St 


c(qt ) - k(X t ), 


Xt — u 0 + 



u s ds, 


while = [ ( J dfdFt is the gain from the self-financing portfolio traded on the futures market. 
A denotes the set of all admissible controls (it, q, 6). More precisely, we will say that a triplet 
(it, q , 9) is an admissible control if: 


• 9 - (<?t)te[o,T] and u = (ut)te\o,T] are adapted processes with values, respectively, in 
[0, q] and [u,u] ; 

• 9 = (#t)tg[o,T] is any predictable real-valued F-integrable process such that the resulting 
wealth is a.s. nonnegative at any time, i.e. 

r 0 + f ir s ds + V t e > 0, t G [0, T\. 

Jo 

The relevant state variable of the problem is Z = ( R , X, D) where R is the wealth of the 
agent, i.e. 

Rt = t’o+ [ Tt s ds + Vf, t G [0, T\. 

Jo 

The dynamics of the state variable is given by 

dR t = [(q t - 1 tt)t/}(D t ) - c(q t ) - k(X t ) + a(t, D t )6 t \ dt + fi(t, D t )9 t dW t , 

dX t = utdt, 

dD t = aD t dt + crdWt■ 

Let us introduce the value function of the optimization problem as 

v(t,r,x,d) = sup E [{Rt) 1 \Zt = (r, x, d)\ 

(u,q,9)eAt 

where At denotes the set of all admissible controls starting at time t. The corresponding 
Hamilton-Jacobi-Bellman equation (hereafter HJB equation) is given by 

— vt — sup Cv = 0, with A := [u,u] x [0, q] x M, (4.6) 

(u,q,0)eA 


with terminal condition 

v(T, r, x , d) = r 7 , 


and where 


(4.7) 



The HJB equation rewrites: 


0 = -v t - adv d + k(x)v r - -a 2 v dd 


1 


— sup \ uv x + [(g — u)'ijj(d) — c{q) + a(t, d)6} v r + d) 2 9 2 v [ 


(■ u,q,6)eA 
Rearranging the terms gives 


0 = -v t - adv d + k{x)v r - -cr 2 v d d 


sup 

(■ u,q,e)eA 


[v x — ijj(d)v r )u — ( c(q ) — qtp(d))v r + a(t, d)6v r + —/3(t, d)~9 2 v. 


(4.8) 


We notice immediately from the above HJB equation that the optimal candidate rule for the 
storage management is 

U t — 'MX(l/j(Dt)Vr>V x ) T U ^-(ll>(Dt)Vr<V x )-> (4-9) 

where, to simplify the notation, we dropped the arguments ( t , Rt, Xt, Dt) from the derivatives 
of the value functions v r and v x . Depending on the ratio r] = v x /i)(d) between the marginal 
utility of one unit of storage and the spot price, it is optimal either to buy and store at 
maximum capacity or to withdraw and sell at maximum capacity. One may have thought 
that this ratio should simply be compared to one. This is not the case as it has to be 
compared to the marginal utility of one unit of wealth, v r . As pointed at the end of Section 
3, we see that the storage control depends on the wealth and hence on the trading activities. 


Furthermore, the heuristic computation above confirms the fact that the optimal control for 
production is to produce until the marginal cost of production equals the spot price, i.e. 
equation (3.2). 


Finally, solving the maximisation problem in (4.8) gives the optimal control for the trading 
portfolio as 


6 


* _ 

t ~ 


a(t,Dt ) v r 
(3(t,D t ) 2 v rr 


(t, Rt, X t , D t ). 


(4.10) 


Since it is likely that v rr is negative because U is concave, one recovers the expected result 
that the agent holds a long position if futures prices are exhibiting a positive trend. Moreover, 
9* is similar to the Sharpe ratio, a tradeoff between the expected trend of the futures prices 
compared to their volatility. 

Remark 4.5. To rigorously solve the optimization problem, one should prove that the Cauchy 
problem (4.6, 4.7) admits a unique solution with the required regularity together with a 
verification theorem. This could be achieved using the techniques developed in, e.g., Pham 
(2002 [32]) for quite a large class of multidimensional stochastic volatility models. However 
adapting such a method to our setting would be very technical and would go far beyond the 
scope of this paper. 


5 Conclusion 

In this paper, we develop a parsimonious model of commodity spot and futures prices, that 
can explain the relation between the spot and the futures prices based on arbitrage arguments. 
The argument does not rely on the storability of the commodity. Our setting shows that 
somehow, the possibility to control the production is like the ability to store, albeit for just 
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an instant, as is the possibility to increase of decrease the amount hold in a storage capacity. 
From this point of view, the financial assets and the futures contracts traded on commodity 
markets are not so different from each others. 

The existence of a risk-neutral measure is shown to be a consequence of the finiteness of the 
value function of an agent maximizing his utility. We also prove that the futures price always 
converges to the spot price, whatever the storability property of the underlying commodity 
is. Finally, we briefly discuss the solution of the trading-production problem faced by the 
agent and we show how the different controls can be separated. In particular, the optimal 
storage policy is impacted by the possibility of trading on the futures market. 
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